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Impact of QED corrections to Higgs decay into four leptons at the LHC 
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At the LHC a precise measurement of the Higgs boson mass (if discovered), at the level of 0.1-1%, will be 
possible through the channel gg ^ H ^ 41 for a wide range of Higgs mass values. To match such an accuracy, 
the systematic effects induced by QED corrections need to be investigated. In the present study the calculation 
of 0{a) and higher order QED corrections is illustrated as well as their impact on the Higgs mass determination, 
once realistic event selection criteria for charged leptons and photons are considered. 
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1. Introduction 

Assuming that the Higgs boson will be discov- 
ered at the LHC, a measurement of its mass, for 
a large range of values, with a relative experi- 
mental precision of 0.1-1% will be possible, by 
combining ATLAS and CMS measurements and 
for an integrated luminosity of 300 fb~^ per ex- 
periment p. To match such an accuracy, the 
impact of QED radiative corrections on the Higgs 
mass determination through the gluon- fusion pro- 
cess gg ^ H ^> Al (4/ = 4e,4^, e+e^^+/i^) 
must be evaluated. To this end, exact 0{a) 
and higher-order QED corrections are calculated, 
as described in Sect. 12 and their effect on the 
Higgs mass extraction in the presence of realis- 
tic event selection criteria is evaluated, as dis- 
cussed in Sect. O Although the electromagnetic 
corrections affect only the final state, the effects 
of the production process have to be considered, 
because the typical event selection applied on lep- 
tons are not Lorentz invariant. The complete 
process pp ^ H ^ 4Z(7) has been simulated in 
the narrow width approximation, well justified for 
low Higgs masses, thus allowing to factorize the 
production and decay processes. The obtained 
preliminary results shown in Sect. |21 indicate that 
QED radiation effects should be carefully consid- 



ered in view of the expected precision at the LHC. 



2. The calculation 

QED radiative corrections are calculated on 
top of exact tree-level matrix elements for the de- 
cay H -^ Z^*'>Z^*^ -^ 41, which, at tree level, con- 



'H+H' 



sists of one diagram for the final state e^e 

and two diagrams for equal flavour leptonic pairs 

(see Fig. Pi . 

Since the tree-level is mediated by electrically 
neutral Z-bosons, QED corrections can be safely 
calculated in a gauge invariant way as a subset of 
the complete electroweak corrections T, namely 
neglecting W^-boson exchange contributions. Two 
complementary approaches have been adopted: 
the former is based on the parton shower tech- 
nique, in the realization of Ref. j3j, allowing to 
calculate QED corrections in the leading loga- 
rithmic approximation both at 0{a) and to all 
orders; the latter relies upon a complete 0{a) 
perturbative calculation. While in the first case 
the correction is completely factorized over the 
tree-level, in the diagrammatic calculation some 
care has to be devoted to the treatment of the 
Z-boson width in the virtual and real corrections 
in order to safely control the infrared (IR) can- 
cellations. As for any one-loop calculation, the 
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expression for the decay width, fully differential 
over the final state, can be written as 



dr = {dr)B + {dr)v + (dr)^. 



(1) 



where the subscripts stand for Born (B), Virtual 
(V) and Real (R) contributions. The IR diver- 
gence is regularized by means of a small photon 
mass A. The real corrections are calculated ana- 
lytically in soft approximation for A < E^ < fco 
(where kg is the soft-hard separator), and by 
means of exact matrix elements for E^ > fco, 
with finite fermion masses and A = 0. The 
real hard photon emission diagrams (one of them 
is shown in Fig. ^ have been calculated ana- 
lytically with FORM 4 and cross-checked with 
ALPHA [S] . The virtual corrections consist of ver- 
tex and self-energy diagrams (see Fig. I^J, which 
are symbolically written in terms of Passarino- 
Veltman form factors [H] and evaluated numeri- 
cally with LoopTools [7|. In addition, also box 
and pentagon diagrams are present. One exam- 
ple of five-point graphs is given in Fig. 0] The 
pentagon diagrams are reduced, with the help 
of FORM :4, to combinations of four-point form 
factors by means of the techniques introduced 
in Ref. _8_, in order to avoid numerical instabil- 
ities due to vanishing Gram determinants. The 
method has already been successfully used for the 
calculation of the 0{a) electroweak corrections to 
e~^e^ — > 4 fermions, where also six-point func- 
tions are involved |9j. Adopting the approxima- 
tion of vanishing fermion masses whenever pos- 
sible in the virtual corrections and performing 
the calculation in the 't Hooft-Feynnian gauge 
^ = 1, the involved five-point functions are at 
most of rank two and therefore are free of ul- 
traviolet divergences. An additional complica- 
tion is due to the presence of the unstable Z 
bosons. In order to avoid singularities in the 
phase space, the introduction of the Z width is 
required, which could spoil the IR cancellation 
between virtual and real corrections. In fact the 
IR divergences, contained in the non-factorizable 
five-point diagrams, are cancelled by the interfer- 
ence between real (tree-level) radiation from dif- 
ferent external legs. A solution is given by the 
"complex mass scheme" (introduced in Ref. [TU] 
for lowest order processes and generalized for one- 





Figure 1. Born diagrams for the decay H -^ 4e. 



loop calculations in Ref. [H]), where the Z mass 
is shifted on the complex plane with fixed width 
A^l -^ Af| - iTzMz, both in tree-level and in 
loop diagrams and the couplings become complex 
quantities, in order to respect the Ward identi- 
ties. Considering that self-energies and vertex 
corrections, neglecting terms of 0(m?/Q^), are 
already factorized over the tree-level, and that 
with complex Mz the IR singularity can be fac- 
torized over the tree-level also for five-point di- 
agrams, the 0{a) QED corrected Higgs partial 
width can be written as: 



dVoia) = (dr)Bx(i + ,sr* + '5r ) + (^r)fl^ 

+ [idr)^^-idT)BxS'y''% (2) 



where Sy^^ refers to the contribution of self- 
energies and vertices, (5y~ refers to the IR scalar 
three-point functions representing the IR part of 
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Figure 2. The real-radiation diagrams are ob- 
tained attaching a real photon leg to every final 
state fermion. 



the five-point functions and {dT)y is the com- 
plete contribution of the pentagon diagrams. By 
construction, the IR divergence has been factor- 
ized over the tree-level, thus allowing for a con- 
sistent IR cancellation with the real part {dr)ii, 
and the remainder [{dr)y - {dT)B x 6y^'^] is free 
of divergence. While the numerical evaluation of 
two- and three-point scalar functions with com- 
plex masses is performed with LoopTools v2.2, 
the four-point scalar scalar functions is reduced 
to unidimensional integrals which are evaluated 
by means of the adaptive integration package 
CUBA ^I]. With the same numerical algorithm 
several checks have been performed of the two- 
and three-point functions with complex masses 
obtained with LoopTools (which uses the formu- 
lae implemented in FF |12j). finding up-to-digit 
agreement. 

In order to estimate realistically the impact of 
QED corrections on the Higgs mass determina- 
tion, the complete production and decay process 
pp -^ H ^ Al + {n-f) has to be simulated, con- 
sidering typical realistic event selection criteria. 
Since for a Higgs mass value below the real Z-pair 
production threshold the total width is very small 
compared to its mass, the narrow width approx- 
imation is adequate. In this approximation the 
calculation can be split in on-shell Higgs produc- 




Figure 3. Vertex (upper panel) and self-energy 
(lower panel) diagrams contribute to the factor- 
izable corrections. 



tion X decay, according to the following formula 



a{pp -^ H ^ 4:1 + {wy)) = 
/ dxidx2fgixi, p)fg{x2, p) 



^gg^H S{xiX2S - MJj) 

dT{H ->4:l + (n-f)) 



9 (cuts). 



(3) 



where fg{x, /j,) are the gluonic parton distribution 
functions, agg^H is the parton level Higgs pro- 
duction cross section and 0(cuts) is a step func- 
tion accounting for the experimental cuts. 

A Monte Carlo code has been developed, based 
on the event generator ALPGEN ^3j for Higgs pro- 
duction and an original library H24F for the decay 
into four leptons with QED radiative corrections, 
taking into account realistic event selections for 
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Figure 4. Example of pentagon diagram, which 
contributes to the non-factorizable corrections. 



charged leptons and photons. The additional ap- 
proximation of neglecting the Higgs transverse 
momentum is assumed in the present study. A 
more detailed investigation will be presented else- 
where. 

3. Numerical results 

In order to quantify the shift induced by QED 
corrections on the Higgs mass determination, 
binned x^ fits to the (four leptons) invariant mass 
distribution (see Fig.EJ have been performed, fol- 
lowing the strategy described in detail in Ref . ^H . 
A reference invariant mass distribution is pro- 
duced at Born level for a given Higgs mass and 
N invariant mass distributions are produced in- 
cluding radiative corrections for different Higgs 
masses, around the reference value (with 20 MeV 
spacing). N x^ values are then calculated, ac- 
cording to the following equation 



X^imn) 
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I— bins 



\ CTQED O-Born J 



A^ 



<^QED 



d(7i 



(4) 



where dai^Bom/crBom and dai^QEo/o-QED are the 
Born-level and QED corrected cross sections rel- 
ative to the i-th bin (normalized to the corre- 
sponding integrated cross sections). The Higgs 
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Figure 5. Higgs invariant mass in the laboratory 
frame for the decay H —f e"*"e~^^/i~, after apply- 
ing experimental cuts and momentum smearing, 
as obtained at tree level, with Parton Shower and 
with exact 0{a) corrections. 



mass shift is derived from the minimum of the x^ 
distribution. 

In order to be close to the experimental selec- 
tion, lepton identification requirements and kine- 
matical cuts, as well as uncertainties in the energy 
and transverse momentum measurements of the 
photons/leptons in the detector, have been imple- 
mented according to Ref. [Q . The results of this 
preliminary analysis are given in Tab. ^ show- 
ing the Higgs mass shifts due to 0{a) (second 
column) and higher-order (third column) correc- 
tions, for the three possible four lepton decays 
and for a Higgs mass of 130 GeV. It can be seen 
that the mass shift due to multiple photon ra- 
diation is, as a rule of thumb, of the order of 
10% of that caused by one photon emission and, 
therefore, non-negligible in view of the expected 
precision at the LHC. Furthermore, it has been 
preliminarily observed that, at the present level 
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Table 1 

The Higgs mass shifts due to 0{a) and higher-order QED corrections to gg 

different lepton final states and for a Higgs mass of 130 GeV. 
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of accuracy of the investigation, exact 0{a) and 
leading logarithmic 0{a) QED corrections induce 
the same mass shifts. A more detailed analysis, 
as well as a presentation of further numerical re- 
sults, is left to a future work. 
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